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^ ABSTRACT 

^ The majority of debris discs discovered so far have only been detected through infrared 

excess emission above stellar photospheres. While disc properties can be inferred from 
unresolved photometry alone under various assumptions for the physical properties 
of dust grains, there is a degeneracy between disc radius and dust temperature that 
depends on the grain size distribution and optical properties. By resolving the disc we 
can measure the actual location of the dust. The launch of Herschel, with an angular 
H resolution superior to previous far-infrared telescopes, allows us to spatially resolve 

CS| more discs and locate the dust directly. Here we present the nine resolved discs around 

A stars between 20 and 40 pc observed by the DEBRIS survey. We use these data 
to investigate the disc radii by fitting narrow ring models to images at 70, 100 and 
160 fim and by fitting blackbodies to full spectral energy distributions. We do this with 
the aim of finding an improved way of estimating disc radii for unresolved systems. 
The ratio between the resolved and blackbody radii varies between 1 and 2.5. This 
ratio is inversely correlated with luminosity and any remaining discrepancies are most 
likely explained by differences to the minimum size of grain in the size distribution 
or differences in composition. We find that three of the systems are well fit by a 
narrow ring, two systems are borderline cases and the other four likely require wider 
or multiple rings to fully explain the observations, reflecting the diversity of planetary 
systems. 

Key words: circumstellar matter - planetary systems - infrared: stars - stars: indi- 
vidual: (a CrB, Uma, A Boo, e Pav, C Eri, 7 Tri, p Vir, 30 Mon, Tri) 
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1 INTRODUCTION 

Debris discs have now been discovered around hundreds of 
stars. In most cases we only know of their existence due 
to observations in the infrared showing a greater flux than 
is expected for the star's photosphere alone, which is at- 



© 0000 RAS 



2 M. Booth et al 



Table 1. Star information. The first column gives the Unbiased Nearby Stars (UNS) identifier from the DEBRIS survey target 
list jPhillips et al.|2010| ). The first letter gives the spectral type of the primary and the number gives the rank order of distance 
from the Sun. 



UNS ID 


HD 


Name 


Spectral Type 


Distance (pc) 


Age (Gyr) 


RA 


DEC 


A018 


139006 


a CrB 


A1IV 


23.0 


0.27±0.1 


15 h 34 m 41?268 


+26 42'52'89 


A024 


95418 


P Uma 


AlIVps (Sr II) 


24.5 


0.31±0.1 


ll h 01 m 50?477 


+56°22 / 56y73 


A053 


125162 


A Boo 


A3Va 


30.4 


0.29±0.1 


14 h 16 m 23!019 


+46°05 / 17y90 


A061 


188228 


e Pav 


AOVa 


32.2 


0.25±0.1 


20 h 00 m 35!556 


-72 54'37{'82 


A064 


20320 


C Eri 


kA4hA9mA9V 


33.6 


0.80±0.1 


03 h 15 m 50!027 


-08°49 / liy02 


A067 


14055 


7 Tri 


AlVnn 


34.4 


0.16±0.3 


02 h 17 m 18!867 


+33 o 50 ; 49y90 


A076 


110411 


p Vir 


A3Va 


36.3 


0.50±0.1 


12 h 41 m 53!057 


+10°14 / 08"25 


A082 


71155 


30 Mon 


AOVa 


37.5 


0.17±0.3 


08 h 25 m 39!632 


-03°54 / 23yi2 


A086 


13161 


P Tri 


A5IV 


38.9 


0.73±0.3 


02 h 09 m 32!627 


+34°59 / 14y27 



Spectral types and distances are from [Phillips et al.| ( |2010| . Ages are from[Vican (2012). Right ascensions and declinations are 
the J2000 epoch ( |van Leeuwen|2007| >. 



tributed to dust grains produced by collisions of planetesi- 
mals. Knowing the radial structure of the parent planetes- 
imals is key to probing many aspects of debris discs and, 
by extension, the planetary system as a whole. When no 
spatial is information available, fitting a blackbody to the 
spectral energy distribution (SED) provides an estimate of 
the dust temperature that can be converted into a distance 
from the star. However, a blackbody is a poor approxima- 
tion for small dust grains, which are inefficient emitters at 
long wavelengths and are generally hotter than blackbodies 
(e.g. Bac kman &; Paresce|19 93). 

So far, only a small fraction of these disc candidates 
have been spatially resolved. For a limited number of sys- 
tems, scattered light from the dusty debris can be seen when 
the light from the star is carefully removed and/or blocked 
by a coronagraph, though the interpretation of these data 
are complicated by the fact that the distribution of the 
submicron scattering grains is strongly affected by radia- 
tion forces. At longer infrared wavelengths, thermal emission 
from the dusty debris can be detected and contrast with the 
star becomes less problematic. Moreover, far-infrared obser- 
vations probe grains that are above the size at which grains 
are blown out of the system by radiation pressure and are 
thus expected to trace the locations of the parent planetes- 
imals. 

Together with the SED, resolved images allows us to di- 
rectly measure the radial location of emitting dust, and pro- 
vide constraints on the grain size and composition. Studies 
based on a combination of observations at different wave- 
lengths and using different methods have suggested that re- 
solved debris disc radii are 1-5 times larger than inferred 
using blackbody fits (e.g. Rodriguez and Zuckerman 2012). 
In addition, resolved images can show us radial and vertical 
asy mmetries (e.g.|Kalas &; Jewitt|1995| , azimuthal structure 



has two key programmes dedicated to searching for debris 
discs: Disc Emission via a Bias-free Reconnaissance in the 



(e.g. |Greaves et al.||2005 ) and inner edges (e.g. Boley et al. 
2012), which in turn can be used to infer the existence of 
planets in the system, as was successfully done with ft Pic 
( Mouillet et al.||1997| |Lagrange et al.||2010| and Fomalhaut 
( Kalas et al.|2005| |Quillen||2006| ). " 

Space based observatories have been at the forefront 
of debris disc observations ever since the first, Vega, was 



Infrared/Submillimetre (DEBRIS [Matthews et al.|2010| ) and 
DUst around NEarby Stars (DUNES |Eiroa et al.|2010^ Her- 
schel provides higher sensitivity and resolution than previ- 
ously available at far-IR wavelengths, and therefore can spa- 
tially resolve many debris discs for the first time. 

In this paper we present images (section [2]) of nine de- 
bris discs around A stars from the DEBRIS survey (a CrB, 
P Uma, A Boo, e Pav, C Eri, 7 Tri, p Vir, 30 Mon, P Tri), 
the cold components of which have never been previously 
resolved. These 9 stars constitute the full sample of resolved 
A star discs between 20 and 40 pc. In section [3] we present 
the spectral energy distributions of these systems. In section 
H] we model the resolved images as narrow rings. In section 
[5] we use our results from a consistent fitting procedure to 
compare the blackbody radii and resolved radii. In section 
[6] we summarise the results for each system individually. In 
section [7] we discuss differences between systems that are 
well described by a narrow ring and those that are radially 
broad. The conclusions are presented in section [8] 



2 OBSERVATIONS 
2.1 PACS 

Herschel's PACS (Photodetecting Array Camera and Spec- 
Poglitsch et al. 2010| instrument observes either 



trometer 

at 100 and 160 /xm simultaneously or at 70 and 160 /xm si- 
multaneously. The DEBRIS survey observed all target stars 
at 100 /im and 160 /xm. PACS 70 /im provides little im- 
provement in terms of sensitivity over MIPS (The Multi- 
band Imaging Photometer for Spitzer, Rieke et al. 2004) 



discovered with the Infrared Astronomy Satellite (Aumann 
et al.|1984 ). The Herschel Space Observatory (Pilbratt et al. 



2010 hereafter Herschel) is the latest infrared satellite and 



70 /im and so was not used for most of the stars, however 
the increased resolution made it worthwhile observing at this 
wavelength for resolved systems. Doing so gave us an extra 
160 /xm image, thus making the 160 pun images of these discs 
twice as deep as the majority of the DEBRIS images. Ob- 
servations were performed in scan-map mode at a scan rate 
of 20" /s for a total observing time of 890s (giving us a total 
observing time of 1780s at 160 pan). Eight, Belong scans are 
performed per observation in two directions that differ by 
40°. 
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a Boo 70|Lim 




P Uma 70|Lim 



X Boo 70|Lim 




8 Pav 70|Lim 



a Boo 100(xm 



a CrB lOO^m 



m 




(3 Uma 100|Lim 



XBoo 100|Lim 





a CrB 160|Lim 
/■V 



P Uma 160|Lim 





Figure 1. PACS observations for the stars presented in this paper. 70/xm images on the left, 100/xm images in the middle and 160/i,m on 
the right. The top row shows the PSF star a Boo. The field of view of each image is 3Q"-y&Q" . N is up and E is left. The lowest contour 
represents 3-cr significance and the contours rise in increments of 12-a for 70/im, 6-cr for 100/xm and 3-cr for 160/xm. 
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Figure [T] (continued) 
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Table 2. Observation data. 



Star 



Instrument 



Obs ID 



Obs Date 



a CrB 


PACS 70/160 
PACS 100/160 


1342223846/7 
1342213794/5 


2011- 
2011- 


-07-09 
-02-07 


P Uma 


PACS 70/160 
PACS 100/160 
SPIRE 


1342220802/3 
1342197015/6 
1342195687 


2011- 
2010- 
2010- 


-04-29 
-05-24 
-04-29 


A Boo 


PACS 70/160 
PACS 100/160 
SPIRE 


1342221869/70 

1342210928/9 

1342222905 


2011- 
2010- 
2011- 


-05-29 
-12-09 
-06-21 


e Pav 


PACS 70/160 
PACS 100/160 
SPIRE 


1342220576/7 
1342204276/7 
1342229226 


2011- 
2010- 
2011- 


-05-04 
-09-10 
-09-22 


C Eri 


PACS 70/160 
PACS 100/160 


1342227302/3 
1342216125/6 


2011- 
2011- 


-08-24 
-03-05 


7 Tri 


PACS 70/160 
PACS 100/160 
SPIRE 


1342237466/7 
1342223876/7 
1342237503 


2012- 
2011- 
2012- 


-01-13 
-07-10 
-01-14 


p Vir 


PACS 70/160 
PACS 100/160 
SPIRE 


1342222452/3 
1342212660/1 
1342222904 


2011- 
2011- 
2011- 


-06-10 
-01-15 
-06-21 


30 Mon 


PACS 70/160 
PACS 100/160 
SPIRE 


1342220306/7 
1342196127/8 
1342195715 


2011- 
2010- 
2010- 


-05-10 
-05-10 
-05-01 


P Tri 


PACS 70/160 
PACS 100/160 
SPIRE 


1342237390/1 
1342223650/1 
1342237504 


2012- 
2011- 
2012- 


-01-12 
-07-04 
-01-14 



The data were reduced using HIPE (Herschel Interac- 
tive Processing Environment, version 7.0, flux calibration 
pixel scale of l"/pix for 70 /am and 



FM6, 



Ott 



jgOlO l to 

100 /xm and 2 / 7pix for 160 /im. Nearest-neighbour map- 
ping has been used throughout this paper rather than driz- 
zle mapping to avoid the correlated error issues present in 
the drizzle method ( |Fruchter &; Hook||2002 |Kennedy et al. 
2012a). 



As the observations are performed in scan-map mode, 
the coverage of each point in the map varies. This is es- 
pecially true of data in the DEBRIS survey where a non- 
standard technique has been developed to allow us to in- 
clude data taken whilst the telescope was turning around. 
Such variation in coverage also means that the noise varies 
greatly over the map. Due to this, the pixel-pixel noise in 
the sky background cannot be measured by measuring the 
variation in pixels in an annulus around the source as this 
will include pixels of different coverage. Instead we separate 
the background pixels into 10 bins of coverage and calculate 
the pixel-pixel noise for each bin. 

The FWHM sizes for PACS at 70, 100 and 160 fim 
are 5.6", 6.8" and 11.4" respectively. To discover which of 
the DEBRIS targets have discs larger than these sizes and 
so can be considered to be resolved, Gaussians can be fit 
to all the systems and the FWHMs are compared to the 
PSF FWHMs. As a secondary check we also fit PSFs to all 



13 
12 

I 11 

o 
E 

E 9 
o 

o 8 



Calibration stars * 
Resolved discs □ 



7 8 9 10 
70/xm major FWHM 



1 1 



12 



Figure 2. Gaussian major FWHM for the calibration stars and 
for the stars with resolved discs modelled in this paper. It is clear 
that these stars are greatly extended compared to the calibration 
stars. 



systems and check for significant residuals after PSF sub- 
traction. It is important to also consider the possibility of 
background confusion in these maps, which is discussed in 
more detail in section [231 

Amongst the 83 A stars in the DEBRIS sample, 7 are 
within 20 pc, one of which is resolved - P Leo - and is 

(12011b. 57 of the stars 



modelled in detail by Churcher et al 



are between 20-40 pc and 9 of these host resolved discs. 
These 9 stars form the basis of this paper. Images of the 
nine systems presented here as well as one of the PSF stars 
(a Boo) are shown in figure [l] These have been cropped to a 
region of 39 // x39 // and centred on the peak of the emission. 
A comparison between the FWHMs of these resolved discs 
and the PSF stars is shown in figure [2] Details of these stars 
are given in table [l] and details of the observations are given 
in tabled] 

The fluxes are calculated through aperture photome- 
try and are given in table [3] PACS data are subject to 
fairly strong low- frequency (1//) noise and so the data were 
pre- filtered to remove this. However, the PACS beam has 
been found to extend to large angular scales with wl0% of 
the energy beyond 1'. Some of the flux is therefore filtered 
out during the high-pass filtering. To account for this, the 
fluxes need to be increased by 16=b5%, 19=b5% and 21±5% 
at 70, 100 and 160 /mi respectively (Kennedy et al. 2012a). 
Although this correction factor has been derived for point 
sources, it should also apply to resolved sources as long as 
the source extent is less than the size scales from which 
flux is removed from the beam due to filtering. A further 
issue with PACS data is that the PACS beam exhibits sig- 
nificant variation at 70 fim as described in more detail in 
Kennedy et al. (in prep). Because of this, we have five dif- 
ferent calibration stars that we can use when comparing to 
our observations. 



2.2 SPIRE 

Those systems where a blackbody fit to the spectrum pre- 
dicted sub-mm fluxes that we expected to be detectable by 
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Table 3. Observed fluxes. F G 5 S (mJy) is the observed flux measured by aperture photometry for resolved images and PSF 
fitting for unresolved images (see the aperture diameter column). F* (mJy) is the flux of the star. F^^ sc (mJy) is the excess flux. 
FWHM ma j is the Full Width Half Maximum of the major axis of a fitted Gaussian. 



otar 


A yfim ) 


robs (.nijyj 




r disc (mJyj 


Aperture 


uiameter j 


rpAAfUA/r . (tt\ 

t wmvi ma j { ) 


a CrB 


70 


515.0±25.2 


105.50±1.73 


409.5±25.3 


30 




6.57±0.03 




100 


235.1±12.6 


51.21±0.84 


183.8±12.6 


30 




7.72±0.04 




160 


67.6±2.4 


19.70±0.32 


47.9±2.5 


PSF fit 




10.83±0.29 


(3 Uma 


70 


393.0±19.4 


93.65±1.58 


299.4±19.5 


30 




6.59±0.03 




100 


189.2±9.6 


45.45±0.77 


143.7±9.6 


30 




7.72±0.06 




160 


58.3±10.5 


17.49±0.30 


40.8±10.5 


30 




10.19±0.21 




250 


18.3±3.9 


7.05±0.12 


11.2±3.9 


PSF fit 








350 


14.1±4.3 


3.55±0.06 


10.5±4.3 


PSF fit 








500 


9.1±4.7 


1.72±0.03 


7.4±4.7 


PSF fit 






A Boo 


70 


345.3±17.3 


22.61±0.40 


322.6±17.3 


40 




8.11±0.05 




100 


272.1±15.4 


11.01±0.19 


261.1±15.4 


40 




9.75±0.07 




160 


142.4±12.1 


4.25±0.07 


138.1±12.1 


40 




12.90±0.19 




250 


50.7±5.1 


1.72±0.03 


48.9±5.1 


PSF fit 








350 


21.3±5.3 


0.87±0.02 


20.5±5.3 


PSF fit 








500 


4.2±4.9 


0.42±0.01 


3.8±4.9 


PSF fit 






e Pav 


70 


63.1±4.8 


19.55±0.29 


43.5±4.8 


30 




7.90±0.17 




100 


41.9±4.6 


9.48±0.14 


32.4±4.6 


30 




9.01±0.26 




160 


23.6±3.2 


3.65±0.05 


19.9±3.2 


PSF fit 




13.06±0.76 




250 


5.0±3.9 


1.47±0.02 


3.5±3.9 


PSF fit 








350 


0.7±4.8 


0.74±0.01 


-0.0±4.8 


PSF fit 








500 


0.0±4.7 


0.36±0.01 


-0.4±4.7 


PSF fit 






C Eri 


70 


93.9±5.8 


17.40±0.43 


76.5±5.8 


30 




7.24±0.13 




100 


84.1±5.9 


8.46±0.21 


75.6±5.9 


30 




10.27±0.24 




160 


42.1±0.8 


3.26±0.08 


38.8±0.8 


30 




13.38±0.52 


7 Tri 


70 


777.6±38.8 


21.20±0.40 


756.4±38.8 


40 




10.24±0.04 




100 


718.8±35.5 


10.29±0.20 


708.5±35.5 


40 




11.63±0.04 




160 


444.3±10.4 


3.96±0.08 


440.4±10.4 


40 




14.09±0.08 




250 


186.6±13.6 


1.60±0.03 


185.0±13.6 


PSF fit 








350 


78.1±7.2 


0.80±0.02 


77.3±7.2 


PSF fit 








500 


21.8±5.1 


0.39±0.01 


21.4±5.1 


PSF fit 






p Vir 


70 


230.1±4.3 


11.25±0.20 


218.8±4.3 


40 




7.91±0.06 




100 


154.2±7.0 


5.48±0.10 


148.7±7.0 


40 




9.31±0.10 




160 


67.3±7.0 


2.12±0.04 


65.2±7.0 


40 




13.62±0.39 




250 


37.9±0.8 


0.86±0.02 


37.0±0.8 


PSF fit 








350 


22.7±0.5 


0.43±0.01 


22.2±0.5 


PSF fit 








500 


20.3±0.4 


0.21±0.00 


20.0±0.4 


PSF fit 






30 Mon 


70 


206.2±10.5 


22.08±0.46 


184.1±10.5 


30 




6.52±0.05 




100 


86.5±5.9 


10.71±0.22 


75.8±5.9 


30 




7.77±0.12 




160 


25.5±0.7 


4.12±0.09 


21.3±0.7 


PSF fit 




11.45±0.76 




250 


7.0±3.4 


1.66±0.03 


5.4±3.4 


PSF fit 








350 


4.2±4.5 


0.84±0.02 


3.4±4.5 


PSF fit 








500 


0.7±4.6 


0.41±0.01 


0.3±4.6 


PSF fit 






(3 Tri 


70 


641.2±31.5 


73.01±1.42 


568.2±31.6 


50 




8.38±0.04 




100 


481.1±23.4 


35.47±0.69 


445.6±23.4 


50 




9.87±0.05 




160 


263.6±0.3 


13.68±0.27 


249.9±0.4 


40 




13.81±0.12 




250 


87.1±7.3 


5.53±0.11 


81.6±7.3 


PSF fit 








350 


34.6±5.6 


2.79±0.05 


31.8±5.6 


PSF fit 








500 


5.1±4.9 


1.35±0.03 


3.7±4.9 


PSF fit 
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the SPIRE instrument ( Griffin et al.|2010 ) were followed up 
with observations at 250, 350 and 500 /im. The larger beam 
sizes of 18. 2", 24.9" and 36.3" at these wavelengths meant 
that the systems were not resolved at these wavelengths, 
however these data are still useful for constraining the SED. 
The fluxes are measured by fitting PSFs to the images and 
are given m table [3] 



2.3 Confusion Issues 

At the depths of the DEBRIS survey a large number of 
background galaxies become visible. Analysis of all DEBRIS 
maps shows that the probability of a background source of 3- 
a significance showing up within 10" of the target is about 



Table 4. Stellar parameters from the photo-spheric fitting. 



36% for the 160 /im maps based on |Berta et al. (2011). 
Therefore, it is likely that some of our discs are in fact con- 
fused with background sources. This does not necessarily 
mean that what was thought to be disc emission is all back- 
ground emission, just that care has to be taken when mod- 
elling the disc as has been shown to be the case for (3 Pic 
( |Regibo et aT1|2012| ) and 61 Vir ( |Wyatt et aT1|2012| ) for ex- 
ample. 

Careful analysis of the SED and images for each disc 
shows that one of the discs in our sample is likely affected 
by background sources. The images of p Vir (figure [lj show 
significant emission to the NE at an RA offset of +10 'and a 
DEC offset of +6". In this case the background source can 
be fit with fluxes of 7.2, 11.3, and 17.7 mJy at 70, 100, and 
160 /xm, which are consistent with a submillimetre galaxy of 
redshift, z « 1 ( Blain et al.|2002 ). The aperture for the pho- 
tometry is increased to include both the disc and background 
emission with the background emission then subtracted from 
the aperture flux. 



3 SPECTRAL ENERGY DISTRIBUTIONS 

Debris discs are usually found by excess IR emission above 
the level expected from the stellar photosphere. Thus, the 
first step is to make a model of the stellar spectral energy dis- 
tribution (SED). To do this we start by collating optical and 
near-IR photometry for all of our stars. Because many of our 
stars are very bright, and therefore saturated in the 2MASS 
catalogue (jCutri et al. 2003 |, we use "heritage" UB V ( |Mer- 
|milliod|2006) and ubvy ( [Hauck fe Mermilliod|1998t photom- 
etry, as well as more recent optical and near-IR photometry 
from Tycho 2 ( |H0g et aT 



2000|), Hipparc os (Perryman & 
ESA|[l997|, AKARI 9 fim |lshihara et aL^OlO) , and IRAS 



12 /im ( Moshir et al.|1990 ). For stars that show excess emis- 
sion at near-IR wavelengths, we do not include AKARI and 
IRAS fluxes in the stellar spectrum fitting. 

The stellar photosphere model is then found by least 
squares minimisation of the GAIA grid of PHOENIX mod- 
els (Brott & Hauschildt 2005). Interstellar reddening is not 
taken into account for the photosphere fitting as it is not 
expected to be important for stars within the Local Bubble 



(Su et al. 2005). This fitting process gives us the effective 
temperature and luminosity of the stars (table These 
parameters are also used to interpolate the mass from the 
tables of Schmidt-Kaler ( 1982 ). Comparison of the predicted 



Star 


Teff (K) 


M* (M Q ) 




a CrB a 


9302±61 


2.7 


59.8±1.0 


(3 Uma 


9291±61 


2.7 


59.7±1.0 


A Boo 


8608±50 


2.3 


17.1±0.3 


e Pav 


9794±61 


3.0 


25.6±0.4 


C Eri 


7516±44 


1.7 


10.9±0.3 


7 Tri 


9127±111 


2.6 


25.2±0.5 


p Vir 


8817±55 


2.4 


13.2±0.2 


30 Mon 


9401±84 


2.8 


34.5±0.7 


(3 Tri a 


8011±42 


1.9 


73.4±1.4 



a These systems are close binaries and so the fitted parameters 
are representative of the combined spectra rather than being real 
properties of the stars. 



in the DEBRIS survey that do not exhibit a detectable ex- 
cess at 24 /im shows that the photospheric predictions have 
a 1-a uncertainty of about 2%. Once modelled, the pho- 
tospheric flux can be subtracted from the observed fluxes, 
yielding the disc spectrum. 

With flux density measurements of an excess at just 
two wavelengths a blackbody can be fit to give the colour 
temperature of the excess spectrum. As many of our stars 
have been observed before, more data than just the Her- 
schel observations can be used to more tightly constrain the 
blackbody fit. Photometry from all available sources was col- 
lected. This includes (where available) IRAS (Faint Source 
Catalog v2; |Moshir et aT]|1990|), ISO (|Habing et al.||200l| ), 



Spitzer IRS ( Lebouteiller et al. 2011), Spitzer MIPS (re- 



reduced data from an updated pipeline), WISE preliminary 
release ( jWright et al.|[2010|, AKARI fllshihara et al.pOlO ) 
and SCUBA ( [Williams fe And7ewsl|2006| ). 

From observations at sub-millimetre wavelengths it is 
known that the spectrum falls off more steeply than black- 
body at wavelengths much larger than the grain size due 
to the absorption and emission inefficiencies. To account for 
this, we use a modified blackbody where the standard black- 
body equation is modified by dividing by a factor X\ as 
given by 



il 



(A/Ao) 



A < Ao 
A > Ao 



(1) 



where Ao and ft are free parameters. In some cases there is 
not enough long wavelength data to constrain the fall-off in 
the SED. In these cases values of Ao = 210 /am and /3 = 1 
are used based on the average values for A star discs ( |Wyatt 
et al.||2007| ). The SED fits for the stars are shown in figure 



3| The results of the SED fits are listed in table \E\ 

The wavelengths observed using Herschel are dominated 
by cold dust. Due to the inclusion of shorter wavelength data 
in the SED, it is important to also consider the possibility 
of emission from a warm disc that may be analogous to the 
asteroid belt but could also be evidence of emission from 
dust at a single radial location at multiple temperatures. 
Three systems show evidence for warm emission based on 
their SEDs: p Vir, 7 Tri and 30 Mon. 



photospheric fluxes with MIPS 24 /im photometry for stars 



Morales et al. (2011) showed that the Spitzer data for 
p Vir is best fit by a two ring system with belts at tempera- 
tures of 237 and 74 K. With our increased wavelength range 
for the SED we find this to be broadly correct, with our best 
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fit giving temperatures of 279±107 and 81±3 K. The disc 
around 7 Tri has previously been reported to have tempera- 
tures ranging from 65 K ( Zuckerman &; Song 2004 ; Williams 



fe Andrews|[2006t to 120 K ( [Chen et al.||2006| ~Clearly the 
dust in this system is at more than one temperature so here 
we choose to fit it with two ring components, although we 
will see later that narrow rings are not appropriate for this 
system. Note that the inner components for these stars are 
poorly constrained and improved mid-IR photometry, spec- 
tra and/or imaging will be required to get a better idea of 
the properties of this warm component. Such mid-IR imag- 
ing allowed |Smith & Wyatt| ( |2010| to show that 30 Mon 
(HD 71155) most likely has two components and the inner 
component has been resolved by |Moerchen et "aL (2010) at 



figure [4] (left). A line-of-sight integrator method is then used 
to create the image ( |Wyatt et al.||1999| . For each disc size 
and orientation the best fitting flux is then found for the 
model usingthe least squares fitting code MPFIT (Mark- 
1 r The expected photospheric flux of the star is 



wardt 2009 



about 2 AU, which our value of 3±4 AU is consistent with. 



4 RESOLVED MODELLING 
4.1 Image preparation 

Herschel has a nominal 1-a pointing error of 2.3 ;/ in scan- 
map mod^J although errors as high as 4 /; or more have been 
reportecQ So, rather than centring the images on where we 
expect the star to be, we centre them (to sub-pixel resolu- 
tion) on the peak of the emission. Although this does have 
the potential to introduce correlated noise as we are regrid- 
ding the data, the difference to the final results is insignif- 
icant and this technique has the benefit of allowing us to 
combine the 160 /im images as noted below. Offsets larger 
than the average pointing error may be indicative of asym- 
metric discs or confusion. The RMS value for the offsets of 
these stars is 1.9" and the largest offset is 3.0" and so we 
do not see evidence that the emission is not centred on the 
star's position for any of our targets. 

Since the 70 and 100 /im data are observed at different 
times and each comes with a simultaneous 160 /im observa- 
tion, the 160 /im images are co-added based on the centre 
location of their respective 70 or 100 /im image. 



4.2 Model creation and fitting 

Since debris disc SEDs are generally well fit by single tem- 
perature blackbodies, we begin by assuming that the debris 
discs in each system can be approximated by a narrow belt 
at a distance, R, from the star with a width equal to 0.1R. 
This width is narrower than for some discs generally con- 
sidered to be narrow, such as the asteroid and Kuiper belts, 
but even using a value as high as OAR gives roughly simi- 
lar results. Assuming all the dust is concentrated to a very 
narrow belt also has the benefit of allowing us to ignore the 
degeneracies between the distribution of dust with radius 
and the fall-off of temperature with radius and allows us to 
keep the number of free parameters low. 

Model fitting is done using a grid-based approach. Us- 
ing the parameters listed above, a torus is set up and filled 
with particles. The model disc is then given an inclination 
to the line-of-sight, /, and a rotation, PA ex t, to set it up as 
it would be seen from the telescope. An example is shown in 



1 http:/ /herschel. esac.esa.int/Docs/Herschel/html/ch02. html 

2 http: / /herschel. esac.esa.int / twiki /bin / view/Public /DpKnownlssues 



added in the central pixel and the resulting image is then 
convolved with the PSF. An example is shown in figure [4] 
(right). Note the expected photospheric flux has to be scaled 
down to account for the filtering of the large scale beam (see 
section [2]) . 

For the majority of the stars a Boo was used as the PSF 
reference star. However, there is some variation at 70 /im 
between the PSF stars as noted in section [2] This variation 
can lead to differences in the quality of a fit, so we use the 
PSF that provides the best fit for a particular star. It was 
found that using 7 Dra as the 70 /im PSF reference star for 
a CrB and f3 Uma resulted in an improved fit (see Kennedy] 
et al.||2012b| ). 



In the cases of 7 Tri, 30 Mon and p Vir, the SEDs indi- 
cate the presence of warm dust (see section [3]). For the SED 
fitting we assumed this warm dust to indicate the presence 
of a narrow belt close to the star analogous to the aster- 
oid belt. We could include this as a fixed flux component 
in the image modelling, however, since doing so involves ex- 
trapolating from an already uncertain mid-IR flux to PACS 
wavelengths, we chose to ignore this component. Doing so 
does not have a significant effect on our results as the inner 
components' fluxes are relatively low at these wavelengths. 

This is all then rerun for different values of R, I and 
PA ext . A grid of around 50,000 models was used for each 
system, the exact number depending on how fine a grid was 
needed to accurately determine the errors on the parame- 
ters. 

For each grid model at each wavelength the reduced % 2 
is found 



2 _ E: = p o((fim(i) 

Ared 



F sim (i))/<r) 



(2) 



where a is the pixel to pixel RMS, Fi m (i) and ^sim(^) are 
the fluxes in pixel i of the image and simulation respectively, 
rip is the number of pixels used for the fitting process and 
riparam is the number of free parameters. The best fitting 
model is then the model with the minimum x?ed- The pixels 
used for this calculation are those within 1.5xFWHM of the 
resolved disc (where this is calculated from a Gaussian fit) 
that also have a flux density greater than the 1-a pixel to 
pixel error. To determine whether or not this best fit model 
is in fact a good fit to the data, i.e. whether or not a single 
narrow ring is a good fit, we define a cut-off value for the 
Xred as the value for which there is only a 5% probability 
that such a model is correct. Any model with a x^ed greater 
than this is considered to be incorrect. 

As mentioned in section [2] p Vir has a background 
source very close to the star. If not taken into account, this 



3 Flux was not included as a grid parameter as it is strongly 
constrained by the images so would require a very fine grid to 
accurately determine its value and it is also much less likely that 
using a Levenberg-Marquardt technique to find the minimum % 2 
might encounter a local minima, which is one of the main reasons 
for using a grid method to find the other parameters. 



© 0000 RAS, MNRAS 000, 000-000 



Resolved DEBRIS Discs 9 




1 10 100 1000 1 10 100 1000 1 10 100 1000 

Wavelength, 



Figure 3. SED fits to each of the stars. The solid line represents the combined star and disc fit, whereas the dashed lines represent 
separate fits to the star and disc. The black squares represent the observed photometry and the green diamonds represent the star 
subtracted photometry. The black and green downward triangles represent upper limits for the observed and star subtracted photometry. 
The blue line represents the IRS spectrum. The green dots represent the star subtracted IRS data. No error bars have been plotted for 
the IRS data for clarity, but the uncertainties do increase greatly at shorter wavelengths. 



Table 5. Parameters from a blackbody fit to the SED. Columns 2 to 7 give the temperature, radius and fractional luminosity for 
the inner belt (subscript i), where this is included in the fit, and for the outer belt (subscript o). Ao and /3 are the parameters for 
the modified part of the blackbody (see section pi. Xred & ves the reduced chi-squared of the fit. 



Star 


T, (K) 


Ri (AU) 


f* (xlO- 5 ) 


To (K) 


Ro (AU) 


fo (xlO- 5 ) 


Ao 


P 


X 2 , 
^red 


a CrB 








121.8±3.4 


40.4±2.3 


1.7±0.2 


74.6±6.7 


1.6±0.2 


0.241 


(3 Uma 








119.1±3.6 


42.2±2.6 


1.3±0.2 


72.1±9.3 


1.3±0.3 


0.664 


A Boo 








91.1±1.1 


38.6±1.0 


5.1±0.3 


162.2±41.9 


1.1±0.6 


0.722 


e Pav 








87.0±3.1 


51.9±3.7 


0.5±0.1 


210.0±0.0 


l.OiO.O 


1.313 


C Eri 








80.1±1.1 


39.7±1.2 


2.0±0.1 


210.0±0.0 


l.OdbO.O 


1.434 


7 Tri 


182.5±38.4 


11.7±4.9 


2.15±1.93 


65.6±2.1 


90.2±5.8 


6.4±0.5 


175.2±21.3 


1.1±0.2 


1.164 


p Vir 


278.6±107.1 


3.6±2.8 


1.48±1.99 


80.7±2.6 


43.1±2.8 


5.0±0.6 


71.8±16.4 


0.7±0.3 


1.124 


30 Mon 


391.4±249.3 


3.0±3.8 


1.73±3.81 


107.8±9.3 


39.2±6.7 


2.9±0.9 


64.7±10.1 


1.5±0.1 


0.753 


P Tri 








86.4±1.0 


89.0±2.3 


3.1±0.2 


162.6±21.2 


1.3±0.3 


0.982 



a In these cases the fall off in emission at long wavelengths is unconstrained by the data and so Ao and (5 are fixed to values believed 
to be typical for A stars (see text). 
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Figure 4. Left: Best-fit model of the e Pav disc at 100 /im. Right: Best-fit model at 100 jum of e Pav combined with the stellar photosphere 
and convolved with the PSF. 



can easily influence the fitted parameters. To stop this back- 
ground source from influencing the model fitting, this region 
is masked for the purposes of calculating the x 2 . 

1-a errors on the parameters are then calculated by find- 
ing the maximum and minimum parameters for models with 
a x 2 within A(% 2 ) of the minimum % 2 5 where A(% 2 ) = 3.53 
when there are 3 parameters (see chapter 15 of Press et al.| 
1992). These parameter uncertainties work on the assump- 



tion that the best-fit model is the correct description of the 
disc. Therefore, they do not take into account disc geome- 
tries different to a narrow ring or the real optical properties 
of the grains and so the quoted uncertainties are likely to 
be somewhat underestimated. For this reason, in the cases 
where a disc has a best-fit model with a x 2 too high for 
the model to be considered believable, we do not provide 
parameter errors as they would be meaningless. 

The parameters for the best-fit models are shown in 
table [6] This shows the best fits to each of the wavelengths 
individually and the grid model that best fits all of the wave- 
lengths is shown in the rows marked 'all'. The residuals left 
after the 'all' best-fit model has been subtracted from the 
observations are shown in figure [5] Discussion of individual 
systems is left for section [6] 

4.3 Problems with a single narrow ring fit 

In this paper we have been focussing on fitting all the sys- 
tems with a single, narrow ring to keep the model simple 
and avoid any need to make assumptions about the grain 
properties. This is seen to work well for spectral energy dis- 
tributions. The same can be said for the case of /3 Uma where 
the fits to different wavelengths result in similar best fitting 
radii. However, in most cases the individual images show an 
increase in best fit radius for increasing wavelengths. Since 
dust further from a star will be colder it will contribute 



more at longer wavelengths. A wide disc seen at low signal- 
to-noise will appear as rings with increasing radius for in- 
creasing wavelength. To test whether this is the case for any 
of our systems we have also calculated which model best 
fits all images simultaneously (table [6]) . The residuals and 
Xred f° r this best fit show a narrow ring to still be consistent 
with the data for e Pav and 30 Mon. This mult i- wavelength 
best fit is dominated by the 70 fim image due to its higher 
signal-to-noise ratio and better resolution. For ( Eri and p 
Vir the 

Xred i s above the cut-off value but the residuals are 
still small and we consider these borderline cases that might 
be affected by an incomplete understanding of the noise. 

Other systems show more obvious deviations from a sin- 
gle narrow ring. Judging by the Xred m table [6] and the resid- 
uals in figure [5] we can see that a CrB, A Boo, and /3 Tri are 
poorly fit at 70 /am and 7 Tri is poorly fit at all wavelengths. 
This implies that these systems have wide belts and/or mul- 
tiple narrow rings. 



5 DEBRIS DISC RADII 

From the previous sections we now have radii calculated 
from fitting a modified blackbody to the SED and radii cal- 
culated from fitting a narrow ring to the images. Figure 
[5] shows the resolved radii plotted against blackbody radii. 
Here we define a factor F that is the ratio of the resolved ra- 
dius to the blackbody radius. Those systems demonstrated 
to be poorly fit by a narrow ring in section [431 are shown in 
red here and in figures [7] and [8] to emphasize that these are 
not narrow ring systems. We still include them as the narrow 
ring fit will give the location that contributes most to the 
flux in the image and therefore the location of the grains that 
contribute most to the photometry at these wavelengths. As 
many of the systems are not well fit by a narrow ring the 
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Figure 5. Residuals after the model that best-fits all wavelengths has been subtracted from the image. 70/im images are on the 
100/im images in the middle and 160/xm on the right. The contours represent ±2 and 3-cr significance. 
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Figure [Fj (continued) 





■ pTri B " 160^^ 



error bars in these plots represent the range of fitted values 
for the fits at individual wavelengths. 

We can see that T has some variation, with values be- 
tween 1 and 2.5, which compares well with the literature 

( |2007t 



e.g. Rodriguez & Zuckerman 2012). Trilling et al 



study debris discs around binary systems and find that, if 
they assume blackbody grains, the dust for some of the sys- 
tems seems to be in unstable circumbinary regions. If instead 
small grains are assumed, many of the unstable systems can 
be considered stable but a similar number of systems that 
were formerly considered stable now have dust in unstable 
regions. Our finding that there is a range of T values with 
some systems surrounded by dust that is close to blackbody 



and others with dust much hotter than blackbody makes 
it possible to explain the findings of Trilling et al. (2007) 
without requiring the dust to be in unstable regions. 

The distribution of T does show an inverse correlation 
with luminosity, as can be seen from figure [7] This luminos- 
ity dependence is due to the the dependence of minimum 
grain size on luminosity as discussed further in the next sec- 
tion. 

We might also expect some of this variation to be due 
to the different ages of the systems. It is not known for sure 
how debris discs become stirred and produce a collisional 
cascade. If this stirring is delayed (e.g. Kenyon & Bromley 



2004 Wyatt 2008) then this could change the quantity of 
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Table 6. Best-fit parameters from modelling of the resolved images. The errors on the parameters are 1-cr errors given by a 
A(x 2 ) cut. In cases where there are no upper or lower bounds, a narrow ring is not a good fit for that system. (See the text 
for further details.) For each model we show the reduced chi-squared, x^ed an( ^ degrees °f freedom and use these to determine 
whether or not this narrow ring model is a good fit to the observations. 



Star 


A (/im) 


Radius 


Radius (AU) Inclination (° 


PA ext (° E 


Fi/ (mJy) 


x 2 d 


d.o.f. 


Good 






(") 




from face on) 


of N) 






fit? 


a CrB 


70 


2.0 


45.6 


79.7 


170.3 


418.2 


2.213 


279 


N 




100 
160 


2.3 
3.4 


+ 2.4 

52.6 Z17 

1-71-7 Q +17.8 

' ' ,Z -24.0 


+0.0 

90.0 ^ U 3 
90 +0 -° 

yU.U _25.o 


170.3 ±2 - .l 

_ +14.7 

170.3 1 5 3^ 


197.0 
53.4 


1.085 
1.010 


329 
113 


Y 
Y 




all 


2.0 


46.1 


79.7 


170.3 




1.566 


729 


N 


(3 Uma 


70 


1.9 


46.4 


82.4 


113.8 


312.3 


1.233 


270 


N 




100 
i fin 

1DU 

all 


2.1 
z.U 
1.9 


52. U _ 3 3 
48.3 +\l 
46.8 ± 0.9 


m + -15.5 
in Q 

oq i +u.y 
oy -- L -30.9 

84-3 I4.4 


108.5 ±^; 7 

00 c +32.3 
o8.5 _ 13 5 

in n 

113.8 I2.7 


163.3 

00.0 

- 


0.870 
u. ( yu 
1.041 


287 

54 

649 


Y 
1 

Y 


A Boo 


70 


2.6 


80.1 


51.3 


42.4 


318.9 


1.776 


375 


N 




100 


3.3 


99.7 


41.8 


42.4 


257.4 


1.296 


469 


N 




160 


4.1 


124.2 tl\t 


41.8 +*-* 8 


_ _ +12.9 
55.0 _]_4.7 


144.6 


0.694 


177 


Y 




all 


2.8 


84.5 


48.4 


44.2 


- 


1.891 


1029 


N 


e Pav 


70 


3.2 


104.4 


49.0 


— 10.3 


51.2 


0.904 


214 


Y 




100 
160 
all 


4.0 
6.9 
3.3 


127.9 Z 13 t 
221 4 + 45 - 8 

-49.1 

106.0 ±l 6 Q 


07 9 +8.8 
°' ,z -20.5 

63 6 +26 4 

_23.4 

34.3 ±|:1 


57 q +25.9 
°'- y -25.9 
„ „ „ _i_9n 7 
89.0 i^; 7 

11-4 ±ii-j 


50.3 
31.1 


1.072 
0.912 
1.073 


205 

80 

507 


Y 
Y 
Y 


C Eri 


70 


2.3 


77.9 


42 4 +7 - 9 

-15.9 


25 0+J5-2 

° -U — 15.0 


77.3 


1.061 


216 


Y 




100 


3.8 


127.7 


39.7 


25.0 


81.1 


1.194 


319 


N 




160 


5.3 


1 77 6 +25.0 
_27.0 


45.0 +£* 


25 + 30 - 


50.6 


0.930 


106 


Y 




all 


2.9 


96.5 


45.0 


25.0 




1.290 


649 


N 


7 Tri 


70 


3.7 


126.5 


80.0 


160.5 


770.8 


10.207 


525 


N 




100 


4.3 


148.3 


89.5 


163.7 


717.4 


3.772 


641 


N 




160 


5.0 


171.7 


89.5 


162.9 


450.6 


1.229 


251 


N 




all 


4.0 


1 Q7 A 


OO.Z 


162.1 




6.935 


1425 


N 


p Vir 


70 


2.4 


87.2 


70.9 ±H 


„„ _ + n 

98.51^'g 


236.7 


1.064 


287 


Y 




100 
160 


3.0 
2.4 


109.4 ± 3 4 5 7 
87.2 ±111 


67.1 

5 5 c +29.5 
OO.O _20.5 


90.8 ±H 

00 , +51.9 
t* 6 - L -23.1 


166.4 
76.4 


1.126 
0.901 


284 
72 


Y 
Y 




all 


2.5 


92.3 


69.6 


96.5 




1.174 


651 


N 


30 Mon 


70 
100 
160 
all 


1.8 
2.2 
5.8 
1.8 


67.9 ±H 
83.7 +11 
216.3 ±\%-l 
68.7 


59.2 +jj-6 
43.8 

46.4 +g;S 
56.7 


164.1 + 3 o 6 
167.7 +\™ 

153.3 _ 93 - 3 

167.7 ±™ 


197.7 

88.0 

32.6 


1.024 
0.894 
0.810 
0.978 


244 
220 
85 
557 


Y 
Y 
Y 
Y 


(3 Tri 


70 


2.6 


101.8 


40.3 


66.7 


557.8 


3.343 


447 


N 




100 


3.3 


129.1 


41.3 


63.1 


442.2 


1.568 


543 


N 




160 


4.1 


161.2 ±\% 


44 5 + 4 - 2 

-4.2 


54 7 +6 -° 
oq: -' -6.0 


246.1 


0.869 


226 


Y 




all 


2.8 


109.0 


41.3 


64.3 




3.212 


1224 


N 



small dust in the system with time, which may in turn af- 
fect the radius ratio. Plotting V ratio against age (figure [8]) 
shows no significant change with age and so the effects, if 
any, of delayed stirring would have to occur within the first 
100 Myr. 

5.1 Theoretical explanations for the differences in 
radii 

The reason for the differences between radii inferred by a 
blackbody fit and the actual location of the dust grains is the 



optical properties of the grains. The absorption efficiency, 
Qabs, approaches for A ^> the size of the grain (figure [9| 
and the same is true of the emission efficiencies. This means 
that the grain cannot emit a blackbody at all wavelengths 
and so its temperature increases in order to reach equilib- 
rium temperature. 

The grain size distributions of debris discs are assumed 
to follow a power law size distribution down to the smallest 
size grains in the disc, which is usually considered to be 
the size at which grains are blown out of the system due 
to radiation pressure (Dhi). Therefore the smallest grains 
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Figure 6. Resolved radii versus blackbody radii. The discs in red 
are those for which a narrow ring seems to be a poor fit to the 
disc (see section [4^3} and so the radii simply represents the peak 
of the emission rather than the location of all of the dust. The 
lines represent what we would expect for different values of Y. 
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Figure 8. Ratio of resolved radii to blackbody radii against age. 
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Figure 9. Absorption efficiency for amorphous silicate grains 
based on the core-mantle model of |Li G reenberg ( |1997| . 



Figure 7. Ratio of resolved radii to blackbody radii against stel- 
lar luminosity. 



dominate the cross-sectional area of the system and thereby 
its temperature. The blowout size depends on the luminosity 
of the star and the different luminosities and temperatures 
of stars of different spectral types also affect the temperature 
of a disc for a given location. To demonstrate how a star's 
spectral type affects the T value at different radii we create 
model discs and fit blackbodies to their SEDs, sampled at 
8 wavelengths between 9 and 160/xm, in a similar fashion to 
the SED fitting in section [3] (see also |Bonsor fe Wyatt|2010). 

To do this we use the compositional model of |Wyatt| 
& Dent (2002), which calculates the temperatures, T, and 



absorption efficiencies, Q a bs, using core-mantle grains (Li & 
Greenbe7g||1997 ) . The temperature of the grains are found 
using an iterative process such that 

T(D, R) = ({Q ahs ) Ti I {Q abs ) T(DtR) ) 1/4 T bh (3) 
where Tbb is the black-body temperature and Schmidt-Kaler 



(1982) relations are used to give the temperature and lumi- 
nosity of the star based on an assumed spectral type. The 
absorption efficiency is found using Mie theory, Rayleigh- 
Gans theory or geometric optics in the appropriate limits 
( |Laor feDraine|1993 ) using optical properties from the com- 
positional model. 



The core-mantle model of Li & Greenberg (1997) as- 



sumes the particles to be formed of a silicate core surrounded 
by an organic refractory mantle where the fraction of the to- 
tal volume that is silicate is given by the parameter q s n . The 
silicate material may be either amorphous or crystalline. The 
particles will have a given porosity (p) denning how much 
of the particles' volume is empty space, and a given fraction 
of water (qn 2 o) denning how much of the empty space is 
filled by ice. The size distribution of grains in the disc is 
assumed to follow a power law as n(D)dD oc D~ a dD down 
to a minimum size, D m i n . For the following we assume, un- 
less otherwise stated, the silicates are amorphous, q s n = 1/3, 
p = 0, qu 2 o — 0, 7?= 100 AU, D m in — Du and a — 3.5. 

Figure [To] shows how the radius ratio varies with radius 
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Figure 10. Ratio of resolved radii to blackbody radii against 
actual radii for simulated discs alongside the discs presented in 
this paper. The colour of the symbols and lines represents the 
luminosity of the star. This shows that although we do see a 
correlation between F and spectral type, some differences still 
remain requiring further explanation. 



for stars of different spectral type. This plot also shows the 
results for the resolved discs presented in this paper. The 
model stars use temperatures and luminosities that are the 
average for their spectral type, but significant variation is 
seen in these values for real stars. The symbols and lines 
are, therefore, colour coded by luminosity to ease compar- 
ison. From this plot we can see that the stars with F « 1 
are expected to have these values due to being early type 
stars and so have a large blowout size of ~15 /zm, thus all 
the grains in the system have temperatures close to black- 
body. However, not all the systems fit so nicely with the 
expectations. Further differences are likely to be caused by 
differences in the size distribution. Although finding accu- 
rate fits to every system is beyond the scope of this paper, 
we can demonstrate how we would expect different param- 
eters to vary the temperature of the disc (see also Gaspar 
et al.||2012 ). 

When modelling debris discs it is normally assumed 
that the discs can be represented by a power law slope of 

minimum size 



a = 3.5 (Dohnanyi 1969) going down to 
given by the blowout radius. This may not necessarily be 
the case. Drag forces can cause a turnover in the size dis- 
tribution (e.g. Wyatt et al. 2011) as can interactions be- 
tween blowout grains and the low end of the size distribution 
(Krivov et al. 2000, 2005), thus changing the minimum size 
grain in the distribution. Figure [TT] shows how the radius 
ratio would vary if the minimum size grain in the size dis- 
tribution was different from the blowout size. Since smaller 
grains are hotter than blackbody, increasing the minimum 
grain size reduces the number of small grains and so reduces 
the radius ratio. The effect is stronger for stars of later spec- 
tral type since these stars have smaller blowout grain sizes 
than the early type stars. We also see that it is possible for 
r to exceed unity for an AO star if it has a minimum grain 
size smaller than the blowout size. 

If a more realistic tensile strength curve is used, the 
power law slope of the grain size distribution is steeper than 



Figure 11. Ratio of resolved radii to blackbody radii against 
actual ratio of minimum grain size in size distribution to the 
blowout size for simulated discs. 




Figure 12. Ratio of resolved radii to blackbody radii against the 
power law slope of the size distribution. 



the Dohnanyi relation of n(D)dD oc D ~ a dD ( DurdafeDer- 
lggTllOTBrien fe Greenberg|2003| ). A number of authors 



mott 



have shown that the distribution in debris discs has an in- 



dex of 3.6 or more (Thebault et al.|2003|| 


Krivov et al.|2005| 


Thebault & Augereau|2007| |L6hne et al.j 


2008| |Miiller et al.| 


2010 


Wyatt et al. 2011). Gaspar et al. 


(2012) show that 



this is a robust result for the outer disc zones probed in 
this paper, with a typical index of 3.65=b0.05. For our sim- 
ulations, figure [l2| shows the radius ratio increasing as the 
power law slope steepens, since a steeper slope means more 
small grains in the disc. 

Figure p^l shows how the composition of the particles af- 
fects T. Large values of F are feasible if the grains have no sil- 
icates and no ice, just organics. This results in a peak F value 
of 8 for a porosity of 0.93 and a decreasing F value for porosi- 
ties beyond this as the grain becomes mostly empty space. 
However, it is expected that this kind of organic refractory 
material is formed by UV photoprocessing of ice mantles 



© 0000 RAS, MNRAS 000, 000-000 



16 M. Booth et al 



10 r 



q S !i = °. q H 2o =0 

q S ii=0, q H2 o= 1 

q sil =1/3, q H20 =0 

q s!l = 1/3, q H20 =1 



stars (Tomkin & Popper 1986). More detailed modelling of 
this system is undertaken by (see Kennedy et al.|2012b ). 



6.2 p Uma 

f3 Uma is resolved at 70 and 100 /im. Of all the systems 
in this paper, f3 Uma has the most consistent fitted radii 
across all wavelengths showing the system to be well fit by 
a roughly edge on ring with a radius of 47 AU. This is in 
agreement with the earlier analysis of this data by Matthews 
et al. (2010). This gives it Y = 1.1 as is expected for a star 



of spectral type Al. 



0.0 



0.2 



0.4 0.6 
Porosity 
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Figure 13. Ratio of resolved radii to blackbody radii against 
porosity for grains of different composition in a ring at a distance 
of 100 AU from a star of spectral type A2. The first four models 
are based on the |Li &~G reenberg (1997) core-mantle model and 
the last shows an ice only grain. Note that for the two models 
with qu 2 — 1 5 ah the 'porous space' is actually filled with ice. 
High values of F are possible for organic grains with no silicates 
and no ice. Moderate changes in T are also expected for grains 
with silicates and high porosity. Icy grains can result in values of 
r less than 1. 



on silicate cores (e.g. |Green berg 1982) and so it is unlikely 
that organic only grains would be present in these systems. 
Nonetheless, we present such models here to demonstrate 
the range of values that are possible. For silicate with or- 
ganic grains it is still possible to moderately increase T by 
increasing the porosity so long as no ice is present. For grains 
which only consist of ice it is possible to get values of F less 
than 1 as ice is more efficient at emitting heat than absorb- 
ing starlight. 

Comparison of figures |TT| 1 1 2 1 and 1 1 3 1 shows the minimum 
grain size to be the dominant disc factor in determining T, 
although unusual grain compositions can potentially have a 
large effect. 



6 SUMMARY OF RESULTS FOR INDIVIDUAL 
SYSTEMS 

6.1 a CrB 

a CrB is resolved at 70 and 100 /im and is seen close to edge 
on with a best fit radius of 46 AU giving it T = 1.1. A single 
ring fits the SED perfectly. However, residuals are clearly 
evident in the 70 /im fit, suggesting that the disc may be 



wider than modelled here. Moerchen et al. (2010) observed 



this star at mid-IR wavelengths. Although they do not find 
the excess to be statistically significant, they do find the 
11 /im image to be resolved suggesting a radius of around 
2 AU. Combined with the PACS images this also adds to the 
likelihood that this system is radially broad and the 11 /im 
observation is seeing the inner part of this wide disc. 

This system is an eclipsing binary and the inclination 
of 80° and the position angle of 170° show that the disc is 
consistent with being aligned with the plane of the binary 



6.3 A Boo 

A Boo has a large extent at all wavelengths. A narrow ring 
is a very poor fit to the 70 /im image. When a simulta- 
neous best fit is attempted, 2-a residuals are also seen in 
the 100 and 160 /im images. This shows that a much wider 
belt is necessary to explain these images, probably extend- 
ing between at least 80 and 125 AU. The constraints on the 
inclination show it to be between 40 and 50°. 

This star is the prototype A Boot is star i.e. it has a 
low abundance in refractory elements. It was proposed by 
( Ma rtmez-Galarza et al.|2009 ) that infrared excesses around 
A Bootis stars might be explained by interactions with a dif- 
fuse interstellar cloud, although they noted that a debris disc 
model is more likely for stars like A Boo that are within the 
local bubble. The symmetry of the Herschel images supports 
the presence of a debris disc rather than heated ISM. 



6.4 e Pav 

e Pav is resolved at 70 and 100 /im. Emission at 160 /im 
appears to be offset to both 70 and 100 /im images, however 
this is only significant at the 2-a level and so we cannot draw 
any firm conclusions from it. The results shown in table [5] 
show the fitted radius increasing with increasing wavelength 
from 105 to 220 AU and the inclinations and position angles 
vary between wavelengths. The increase in fitted radii could 
be indicative of a wide disc, however, we find that a narrow 
ring still gives a good simultaneous fit to all the data with 
a radius of 107 AU and inclination of 34°. This results in a 
radii ratio of T = 2.1, somewhat high given its luminosity, 
perhaps suggesting a highly porous or non-silicate compo- 
sition or that some process has led to an increase in grains 
smaller than the blowout size that have yet to be blown out 
by radiation pressure. 



6.5 C Eri 

( Eri is resolved at 70 and 100 /im and marginally at 160 /im. 
The fitted radius is shown to increase from 80 to 175 AU. 
A combined fit at 96 AU with an inclination of 46° results 
in only minimal residuals but the high Xred suggests that 
some extension or a second ring is required to explain this 
disc. The SED is also suggestive of more dust close to the 
star. This star has the highest radius ratio of all the stars 
presented here with T = 2.4, although this is unsurprising 
as it is also the least luminous star in this sample. 
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6.6 7 Tri 

7 Tri has a large extent at all wavelengths and a narrow 
ring gives a poor fit to all images. Multiple temperatures 
are also required to fit to the SED. The fits suggest that the 
disc extends from 125 AU to at least 170 AU and probably 
much further with the disc being close to edge on. 



6.7 p Vir 

p Vir is resolved at 70 and 100 pm with a background source 
present to the NE consistent with a galaxy at redshift z ~ 
1. The fitted radii are consistent at 70 and 160 pm, but 
the 100 pm image requires a larger radius to explain the 
observations and some small but significant residuals are 
present when a combined fit is attempted. The combined 
best fit has a radius of 92 AU and inclination of 20°. This 
gives the disc r = 2.2 and if the dust in the inner disc is 
similar then this should have T ^ 1.5 and an actual radius 
around 4 AU. 



6.8 30 Mon 

30 Mon is resolved at 70 pm and marginally at 100 pm. The 
weak constraints at 100 and 160 pm mean that this system 
is easily fit by a narrow ring with a radius of 68 AU and 
inclination of 56° . This system also has a high mid-IR excess 
with an inner ring expected at around 1 AU ( |Smith &; Wyatt 
|2010| |Moerchen et al.||2010| ). The radius ratio of T = 1.8 is 
higher than expected for a star of this spectral type, perhaps 
indicating a highly porous or non-silicate composition or a 
lower limit to the size distribution than would be expected 
if it were defined by the size of blowout grains. 



6.9 Tri 

/3 Tri has a large extent at all wavelengths. A narrow ring 
is a very poor fit to the 70 pm image and some residuals 
are seen at 100 pm. When a simultaneous best fit is at- 
tempted, significant residuals are seen at all wavelengths. 
Like A Boo and 7 Tri, this disc must be much wider than we 
have modelled here. This system is a spectroscopic binary 
and the inclination of 41° and the position angle of 64° show 
that the disc is consistent with the plane of the binary stars 
dPourbaix||2000|). A more thorough treatment of this debris 
disc and the effect of the binary system on the disc is given 
bylKennedy et al.M2012b). 



7 RADIAL EXTENSION OF DEBRIS DISCS 



As discussed in section |4.3| the images for at least four of 
the systems presented in this paper are indicative of wide 
belts. Kalas et al. (2006) noted a similar distinction between 



narrow and wide belts with some systems having distinct 
outer edges and others with sensitivity-limited outer edges 
and disc widths > 50 AU. 

In our own Solar System the asteroid belt and Kuiper 
belt are both mostly constrained to narrow rings. The main 
asteroid belt is between 2.1 and 3.3 AU and is constrained 
on either side by the presence of Mars and Jupiter. The in- 
ner edge of the Kuiper belt is similarly constrained by the 



presence of Neptune. The outer edge is less well defined with 
high eccentricity objects in the hot classical belt and scat- 
tered disc. However the majority of objects reside between 



^37 a nd 50 AU (based on the CFEPS L7 data, |Petit et al. 
2011) 



The distribution of dust is likely to be wider than this 
due to Poynting- Robertson drag causing dust to drift in to- 
wards the Sun and radiation pressure increasing the eccen- 
tricity of the dust (e.g. |Vitense et aL||2010 ). PR-drag is un- 
likely to be important for the systems presented in this paper 
as they are much more dense than the Kuiper belt and so 
the collision timescale is shorter than the PR-drag timescale 



(Wyatt 2005). Radiation pressure may be more important 



and has two main effects on the distribution. Dust that is be- 
ing blown out of the system forms a blowout grain halo and 
dust that is above the blowout limit is still weakly bound, 
forming a broad disc of eccentric grains. These effects are ob- 
served at these wavelengths around some stars, such as Vega 
flSu et al.|[2Q05j |Sibthorpe et al.||2010| ), but their contribu- 
tion to the flux is small at the longest wavelengths. Blowout 
grain halos cannot be used to explain all wide discs since 
there are cases of wide discs being observed even at millime- 



tre wavelengths (Hughes et al. 



2011 



Wilner et al. 



2012), 



which probe grains much larger than the blowout size. 

There are a number of plausible reasons for distinct 
outer edges in debris discs. For the Solar System, Lykawka 
&; Mukai (2008) have suggested that a planet outside the 



Kuiper belt could be responsible for the outer edge cut-off. 
A similar explanation has been suggested for the ring of 
Fomalhaut (Boley et al. 2012). An alternative dynamical 



explanation could be that objects forming outside this region 
were disrupted by interactions with other stars either in the 
birth cluster ([Adams 2010) or since then ( |Ida et a"T]|2000| . 

Alternatively the edge of a debris disc may be a relic 
of the properties of the original protoplanetary disc. Stellar 
encounters during the early stages of disc formation could 



have stripped material from beyond this distance ( Bate et al 
|2003|) or radiation from a nearby massive star could have 



truncated the disc as is seen to occur in Orion (|Johnstone 
eTaLl[T998l |Mann fe Williams (2010 ) 



Wide discs may simply be discs that formed in the same 
way as narrow rings but were then subjected to dynamical 
interactions with the planetary system. They may be the 
result of dynamical scattering of objects onto eccentric orbits 
similar to the scattered disc component of the Kuiper belt 
but much more massive. However, such events are thought 
to be rare in mature systems (|Booth et al. 2009). Even if 



such an event has occurred, we would have to be seeing 
these discs during or soon after the scattering event as such 
an event rapidly clears material out of a planetesimal disc 
( |Booth et aLp009l ). 

Combining this evidence suggests that these wide discs 
are most likely to be the result of stars that formed with 
large protoplanetary discs that were never truncated. This 
would then require the formation of planetesimals at dis- 
tances up to w 200 AU from the star, further complicating 
planetesimal formation theories that have trouble even ex- 
plaining the formation of Kuiper belt objects. For instance, 
Weidenschilling (2003) suggests that even if the Solar Sys- 
tem's protoplanetary disc was ^100 AU in radius, planetes- 
imals are still unlikely to form beyond 50 AU. Nonetheless, 
recent studies of the formation of planetesimals through 
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gravitational collapse could find ways that allow such ob- 
jects to form at large distances ( Johansen et al.|20b 7, 2012). 



8 CONCLUSIONS 

The increased resolution of Herschel compared to previous 
infrared spacecraft has proven to be a great success in resolv- 
ing debris discs. Observations with Herschel also allow us to 
fill in the spectrum between the mid-IR and sub-millimetre 
i.e. the region where the spectral energy distributions of cold 
discs peak. In this paper we have presented observations of 
nine debris discs observed as part of the DEBRIS survey, 
the cold discs of which have never previously been resolved. 
We fit the images with a narrow ring model to gain an in- 
sight into the parameters of the discs. The discs demonstrate 
the varied nature of planetary systems, ranging from around 
50 AU to 150 AU in size. 

Comparing resolved radii with the radii predicted by 
blackbody fits to the SED shows the resolved radii to be 
larger than the blackbody radii by factors between 1 and 
2.5. Some of this variation can be attributed to differences in 
spectral types of the stars and location of the dust, with the 
ratio of resolved radii and blackbody radii being inversely 
proportional to the luminosity. Remaining differences from 
the expected ratio are most likely to be due to highly porous 
or non-silicate grains or differences to the expected minimum 
size grain in the size distribution, although further work 
is needed to narrow down the exact properties of the size 
distribution. 

We also find that, although a single temperature black- 
body is often a good fit to the SED, a narrow ring is some- 
times not such a good fit to the images. f3 Uma, e Pav and 
30 Mon all appear to be well fit by a narrow ring. £ Eri and 
p Vir show indications that they are probably broader than 
modelled here, a CrB, 7 Tri, A Boo and ft Tri are clear ex- 
amples of discs that cannot be fit by a narrow ring model, 
with very significant residuals left after attempting to fit a 
narrow ring. They are likely to require a broad and/or multi- 
component model similar to 7 Do r (|Broekhoven-Fiene et al.| 
|2012| ) and Leo ( |Churcher etaL|2011) . 

The Herschel Space Observatory has provided us with 
an excellent opportunity to produce a consistent data set 
of resolved discs. Application of the analysis presented in 
this paper to the rest of Herschel's resolved discs will allow 
further exploration of debris disc properties not previously 
accessible with only photometric data and give us further 
insight into how disc properties vary around stars of different 
spectral types. 
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